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ABSTRACT 

At redshifts around 0.1 the CFHT Legacy Survey Deep fields contain some 6 x lO'' galaxies spanning 
the mass range from 10^ to 10^^ Mq. We measure the stellar mass dependence of the two point 
correlation using angular measurements to largely bypass the errors, approximately 0.02 in the median, 
of the photometric redshifts. Inverting the power-law fits with Limber's equation we find that the auto- 
correlation length increases from a very low 0.4/i~^ Mpc at 10^'^ Mq to the conventional 4.5/i~^ Mpc 
at 10^°'^ Mq. The power law fit to the correlation function has a slope which increases from 7 ~ 1.6 at 
high mass to 7 ~ 2.3 at low mass. The spatial cross-correlation of dwarf galaxies with more massive 
galaxies shows fairly similar trends, with a steeper radial dependence at low mass than predicted in 
numerical simulations of sub-halos within galaxy halos. To examine the issue of missing satellites 
we combine the cross-correlation measurements with our estimates of the low mass galaxy number 
density. We find on the average there are 60 ± 20 dwarfs in sub-halos with M (total) > 10^ M© for 
a typical Local Group Af(total)/M(stars) — 30, corresponding to MjLy ~ 100 for a galaxy with 
no recent star formation. The number of dwarfs per galaxy is about a factor of two larger than 
currently found for the Milky Way. Nevertheless, the average dwarf counts are about a factor of 
30 below LCDM simulation results. The divergence from LCDM predictions is one of slope of the 
relation, approximately dN/d\nM ~ —0.5 rather than the predicted —0.9, not sudden onset at some 
characteristic scale. The dwarf galaxy star formation rates span the range from passive to bursting, 
which suggests that there are few completely dark halos. 

Subject headings: Local Group; galaxies: mass function; galaxies: clusters; dark matter 



1. INTRODUCTION 

The detailed physics of whatever constitutes the Cold 
Dark Matter (CDM) may manifest itself in the sub- 
megaparsec regime of the LCDM density perturbation 
spectrum from which low mass dark matter halos orig- 
inate. Low stellar mass galaxies are expected to in- 
habit proportionally low mass dark halos. The avail- 
able evidence supports this view, although dwarf galax- 
ies generally have a higher total mass to stellar light ra- 
tio than higher mass galaxies (reviewed in Mateo 1998). 
Clustering theory predicts that the two-point correla- 
tion of dark matter halos relative to the entire density 
field, that is, bias, decreases for lower masses to become 
slightly less clustered than the mean mass distribution 
(|Kaiser 1984 ISheth fc Torm en 1999') . The easiest clus- 
tering measurement to make, since it has the largest sig- 
nal, is to probe nonlinear clustering where the predic- 
tion is not as clear, since it requires numerical simula- 
tions on scales which are not currently well resolved for 
representatively large cosmological volumes. A special- 
ized case of nonlinear clustering is the distribution of 
dwarf galaxies within the halos of larger galaxies, such 
as the Milky- Way, which is well resolved in n-body sim- 
ulations. As a consequence of the high central density 
of low mass halos, their central regions are very re- 
sistant to dissolution when they merge into high mass 
halos, although their outer regions are tidally stripped 
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away. The outcome is a robust prediction that there 
should be hundreds of low mass sub-halos within the 
halo of large galaxy, which is at le ast a factor of ten 
higher than the vis i ble dwarf galaxies (|Moore et al. 19991 : 
iKlvpin et al. 19991 : ISimon fc Geha 2007t ) of the Milkv- 
Way. 

The problem of relating light to mass for dwarf galaxies 
clouds the comparison to dark matter only simulations. 
The first issue is that dwarf galaxies have been studied 
in detail for a single galaxy, the Milky Way. The cur- 
rent census finds that there are about two dozen known 
Milky Way dwarf galaxies, with many of them having 
been discovered only in the last decade. Under the as- 
sumption of a uniform sky distribution (up to a factor 
of five correction) the total numbers are estimated to be 
about 30 to 50 to My of about - 8 and -4, respectively 
(|Mateo 1998HSimon fc Geha 20071) . 

Measurements of velocities for hundreds of individual 
stars within a goodly number of Milky- Way dwarfs al- 
low an accurate estimate of the mass contained within 
the orbits of the stars for those galaxies. Unfortu- 
nately the mass outside the distribution of the stars 
is essentially unconstrained, even if one makes spe- 
cific assumptions of t he properties of the dark halo 
(jPefiarrubia et al. 20081 ). One surprising outcome is that 
the mass enclosed wit hin the stellar distribution is usu- 
ally close to lO"^ Mq (jMateo 19981: IStrigari et al. 20081) 
over a wide range of dwarf galaxy luminosities. The 
kinematic modelling finds that the central density dis- 
tribution appears to have a core with a universal max i- 
mum value for the central density (|Gilmore et al. 2007f) . 
whereas CDM halos should have cusped central den- 
sity distributions. A dynamical upper limit to the to- 
tal bound mass is that if it becomes too large then their 
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gravitational heating of the thin disk of the Milky Way 
would become unacceptably large. Some authors have 
view ed substructure as a desirable source of disk heat- 
ing (jBenson et al. 200^ . In general, the limit on disk 
requires that if the entire predicted satellite population 
is present, then its radial distribution must be less cen- 
trally concentrated than the over all mass distribution 
()Font et al. 200HlArdi et al. 20031 ). In spite of these un- 
certainties the conclusion remains that the Milky Way 
dwarf galaxies at a total mass of 10^ Mq are a factor of 
ten le ss numerous than t he predicted numbers of sub- 
halos (jStrigari et al. 2007f l. 

Here we use the CFHT Legacy Survey deep imaging 
survey to study a large sample of field dwarfs to comple- 
ment the Local Group studies. Only photometric infor- 
mation is available for each of these small, faint galaxies. 
However, we show that the photometric redshifts turn 
out to be sufficiently accurate that we can reliably place 
these objects into wide redshift bins to measure their 
volume density and their two-point correlation function. 
We show that we can recover standard clustering esti- 
mates for massive galaxies to validate our measurement 
techniques. Then we turn to the cross-correlation of 
low and high mass galaxies. To examine the missing 
satellite problem we combine the radial cross-correlation 
functions with our estimate of the volume density of low 
mass galaxies to derive the numbers of satellites. The 
measurements in this paper use a Hubble constant of 
100/ikms^^ Mpc^^ in a flat ^Im = 0.3 universe. For 
comparison to a specific simulation we use h — 0.73 and 
— 0.24, although these parameter changes do not 
make much of a difference at the low redshifts of this 
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2. DATA AND METHODS 



The deep fields of the CFHT Legacy Survey provide 
[uM,gM,i'M,iM,ZM] images to a depth of im — 26.5 
mag in the Megaprime filter system (similar to the SDSS 
filters). For this study we use the Dl, D2 and D4 fields, 
simply because we wanted an odd number of fields to al- 
low median estimators to be used. The stacking combines 
better seeing exposures leading to stacks having images 
of 0.75" in rM,iM, 0.84" in gM and 0.93" in Uj\/. 
The iM filter is the deepest, with the um and zm images 
having about 2/3 the number of objects that the im cat- 
alogue has. The galaxy counts between im — 26 and 26.5 
mag are rising approximately as log = 0.16iM, slightly 
shallower than the 0. 21 seen in the Hubble Deep Field 
() Williams et al. 19961 ). indicating that incompleteness is 
just setting in. 

The PEGASE/ZPEG spectral energy distribution 
modelling program is used to derive indicative star 
formation rates (SFRs), stellar masses and pho- 
tometric redshif ts for all of the galaxies from 
the photometry (iLe Bo rgne fc Rocca-Vo lmerange 20021 : 
iLe Borgne et al. 2004l: TSullivan et al. 200l). Our mod- 
elling uses the Rana & Basu (1992) initial mass func- 
tion which gave slightly better photometric redshifts than 
alternative IMFs. The particular modelling parame- 
ters adopted limits the maximum specific SFR to about 
10~*yr~^, although this has no effect on our results. The 
modelled stellar masses depend on the IMF chosen with 
alternative IMFs giving systematic mass differences of 
about a factor of two. For instance, the Baldry & Glaze- 
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The redshift distribution of the 10^ - 10'' M© (red) 
10^^ M0 galaxies. The distributions are normalized 



Fig. 1. 

and 10^'' 
Y,n{z)Az = 1. 

brook (2003) IMF in ZPEG gives median galaxy masses 
that are on the average 40% smaller than our adopted 
IMF. All of our clustering and density analysis uses mass 
bins a full decade wide to give good statistics and mini- 
mize the effect of the modelling uncertainties in shifting 
galaxies from one bin to another. 

The Um calibration has been updated from the one 
in Sullivan et al (2006) which helps improve the pho- 
tometric redshift accuracy. Comparison of the photo- 
metric redshifts to spectroscopic redshifts demonstrated 
a redshift ac curacy of a median error of only 0.02, or 
0.012/(l-hz) (|Sullivan et al. 20061) . The redshifts are not 
sensitive to the IMF choice. Redshifts derived from su- 
pernovae in some of the host galaxies provide redshifts 
for galaxies as faint as im — 27 mag. There are no spec- 
troscopic redshifts for the faintest galaxies themselves, 
however the errors should be comparable to the brighter 
galaxies, although almost certainly somewhat larger due 
to increased photometric noise with possible systematic 
errors as a result of the lower metallicities of dwarf galax- 
ies. To minimize the impact of redshift errors on our clus- 
tering measurements we make angular clustering mea- 
surements over redshift bins significantly wider than our 
redshift error. Comparison of the redshifts from different 
parameter choices for ZPEG finds that at least 80% of 
the redshifts, and often more than 90%, will be within 
0.02 of the adopted redshift, with the other distributed 
in long, low amplitude tails. 

The magnitude limit of iab — 26.5 at redshift 0.15 
corresponds to galaxies with absolute magnitudes of 
Mi < —12. Figure [1] shows the redshift distribution of 
galaxies in the range 10^ — 10^ Mq and the rising redshift 
distribution of the more massive galaxies 10^" — 10^^ Mq. 
The low mass galaxy numbers initially rise as would be 
expected for a complete sample, however they begin to 
fall below a uniform density rise of z^Az beyond red- 
shift 0.05 as galaxies with low star formation rates, hence 
lower luminosity at fixed mass, fall out of the sample. 
For a cross-correlation of high and low mass galaxies we 
need to select a redshift range which balances the two 
populations. Given the drop in low mass galaxy num- 
bers beyond z = 0.18 we use that as our upper red- 
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Fig. 2. — A section of the Dl field with 3" circles around 10^ — 
f 0® Mq and 9" circles around galaxies more massive than 10^ Mq 
in the z = O.f — 0.18 range. There are about 25 seeing disks of 
0.75" radius per detected galaxy (at any redshift) in this field. 
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Fig. 3.— The sky locations in field Dl of the 10^'^ Mq and 
"l^glO— 12 galaxies selected to be in the redshift range 0.1 to 
0.18. 



shift. A small section of the sky in the Dl field is shown 
in Figure [2] with galaxies in the z = 0.1 — 0.18 range 
with masses 10^ - 10*^ Mq and 10^ - lO^^ marked. 
In total there are approximately 5 x 10^ detectable ob- 
jects in this field of 0.85 square degrees after masking, 
putting the low mass ga laxies at the confusion limit 
(jScheuer 19571 : iHogg 200l"h . 

The relative positions in the entire Dl field of galaxies 
in the mass ranges 10^ - 10^ Mq and 10^° - 10^^ Mq in 
the selected redshift range are shown in Figure [31 At 
redshift 0.15 the field has a transverse width of about 
10/i^^Mpc. The holes mask out the bright stars along 
with bogus objects and varying background light near 
them. The clustering pattern of the high mass galaxies 
is clearly evident. 

We will measure the clustering with the angular two- 
point cross-correlation function. The two-point func- 
tion gives the sky density of galaxies of type 2 around 
those of type 1 at an angular separation 6 as ^12(6') = 
n2[l + wi2{0)] (the expression is symmetric in the types). 
We use the Landy-Szalay (1993) estimator as modified 
to have the requir ed symmetry in the two populations 
(|Blake et al. 20061 ). For two samples of galaxies, that 
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Fig. 4. — The angular auto-correlation as a function of mass in 
bins of 10^- 10*^, 10^ -lO'', lO'^-lO*, 10** -10^ and 10^ -lO^^ Mq 
(triangles, diamonds, pentagons, hexagons and open circles, or, 
red, green, blue, cyan and magenta, respectively). The errors of 
the mean are displayed. 



share a common distribution area, 



w{9,9 + A9) = 



D1D2 - DiR ~ D2R + RR 
RR ' 



(1) 



where D1D2 is the number of pairs between samples one 
and two that have angular separations between and 
9 + A6'. For the auto-correlation the pairs are measured 
within the single sample of interest. The random sam- 
ple, R, is a uniformly distributed sample of about 10^ 
points, generated with the same boundaries and masked 
regions as the galaxy samples, which in our case have 
the same underlying smooth distribution. The DiR and 
D2R terms are the pair counts with the random points 
and, RR, the random points with themselves. For the rel- 
atively weak clustering visible being measured the clas- 
sical estimator, DD/DR — 1, is substantially less accu- 
rate. Clustering measurements do not depend on sample 
completeness to be accurate, provided that any incom- 
pleteness is not selective in location. 

3. CORRELATION RESULTS 

Figure [H shows the angular auto-correlation of the 
galaxies as a function of stellar mass. We display the 
median (the means are not significant different) of the 
Dl, D2 and D4 field two-point function. The error bars 
are the errors in the means of the three fields. The pair 
counts in a single field for the smallest count bins are 
10^ — 10'^, depending on mass, so the shot noise is usu- 
ally under 10%. The volume of each field is only about 
10*h~^ Mpc^ and the field has a linear span of about 
6h~^ Mpc. The fractional amplitude of the linear CDM 
density perturbation spectrum is such that a sphere of 
8h~^ Mpc radius has a fractional variance of 0.7-0.9, and 
the nonlinear variance will be even larger. Accordingly 
we expect the field-to-field variance in excess of unity to 
overwhelm the relatively small shot noise in our measure- 
ment. 

We fit the power law, w{9) — Au,9^ , to the angular 
clustering measurements. The resulting fits, over the 
angular range 3-400", along with their linear error con- 
tours are shown in Figure \5\ The point to point errors 
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Fig. 5. — The power law fits to the angular auto-correlation func- 
tion. The circle size increases with the mass range of the galaxies 
for mass bins of 10^ - 10*^ (red), 10'^ - 10^ (green), lO'^ - 10** (blue), 
10* - lO'' (cyan) and 10^ - 10^^ Mq (magenta) for the filled con- 
tours. The sizes of the circles centered on each point increases with 
the mass of the objects. 

are sufficiently large that the power law is an accept- 
able fit to the data. Limber's equation inverts the angu- 
lar clustering into a 3D clustering, ^(r) = (ro/r)'^, with 
Aw = J ri?{z)rQr{zY~^ cH{z) dz /[j n{z)dzY, where we 
assume that the clustering is fixed in co-moving co- 
ordinates, r{z). We take the true n{z) to be a con- 
volution of the photometric n{z) with a Gaussian red- 
shift distribution having a az of 0.05 or 0.10, mainly to 
demonstrate that our results are not sensitive to red- 
shift errors. Results for both cr^ are shown in Fig- 
ure [HI We adopt az = 0.05 for our quoted results. The 
error matrix is calculated as the linear derivatives of 
the [Aw , S\ errors into [tq , 7] , which in the log tq axis 
of Fig. [S] become slightly curved ellipses. For large 
masses, M > 10^ M0, the auto-correlation length ro 
and power law index, 7, measured here are in good 
agreement wit h the values measured in field galaxy red- 
shift surveys (iMaddox et al. 1990t iNorberg et al. 20021 : 
IZehavi et al. 2005D . giving us confidence in that the pho- 
tometric redshifts are acceptably accurate. 

The low mass galaxy correlation length is impressively 
low, approximately OAh~^ Mpc with 7 ~ 2 at a median 
redshift of 0.14. If we take the dark matter correlation 
to have an equivalent length of 5h~^ Mpc with the same 
7 then the equivalent bias value, b{M), is 0.1 although it 
is important to recognize that the measurement applies 
only to the non-linear regime. Current simulations do 
not reach such low masses in a field volume. The ex- 
trapolation of the Sheth & Tormen (1999) bias equation 
(their Eq. 12) to low mass at low redshifts corresponds 
to ly approaching zero, or, & ~ 1 -I- (2p — l)/5 — 1, where 
p — 0.3 and 5i ^ 1.68, giving b ~ 0.76, far above what 
we measure for our low mass galaxies. The Sheth & Tor- 
men result was derived for higher masses and for large 
scale clustering of dark halos alone, and are not directly 
testable with our measurements. Hamana et al. (2001) 
address the radial dependence of bias and do find that the 
nonlinear bias for 10^^ M0 halos on a scale of lh~^ Mpc 
does drop into the range of about 0.1 that we measure. 
A clear test awaits better resolved simulations. 
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Fig. 6. — The co-moving correlation length rg and the 3D power 
law, 7, derived using Limber's equation. The photometric redshift 
distribution is convolved with Gaussian redshift errors of 0.05 and 
0.10 (believed to be an over-estimate) in redshift, with the larger 
error always giving a larger ro. The z = 0.1 — 0.18 sample is shown 
as filled contours, the z = 0.1 — 0.35 with an open contour and the 
two bins of 2 = — 0.05 hatched. The mass bins are 10^ — 10^ 
(red), lO'' - 10^ (green), 10^ - 10** (blue), 10* - lO'' (cyan) and 
10^ - 10^2 (magenta). The circle size increases with the mass 
range of the galaxies for mass bins of 10^ — lO'', 10^ — 10^, lO'' — 10*, 
10* - 10^ and 10^ - lO" Mq. 
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Fig. 7. — The angular cross-correlation function of 10^" — 
10^2 M0 galaxies with galaxies in the 10^^ - 10*^ (red), 10^ - 10^ 
(green), lO'^ - 10* (blue), 10* - lO'' Mq (cyan) and 10« - lO"' Mq 
(magenta) range. 

The particular problem of the clustering of low mass 
galaxies in the halo of the Milky Way has been ex- 
amined in consid erable depth as a possible discrepancy 
of LCDM theorv (iKlvpin et al. 1999: 'Moore et al. 19991: 
iDiemand Kuhlen fc Mad au 2007; Springel et al. 2008). 
The simulations find hundreds of sub-halos orbiting 
within a large galaxy size halo which is a factor of at 3 
to 10 more than the number of dwarf galaxies. Our sam- 
ple allows us to investigate the missing satellite problem 
beyond the Milky Way. The first step is to measure the 
cross-correlation between low and high mass galaxies. 

The cross-correlation between two galaxy samples is 
measured with the same methods as we used for the 
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Fig. 8. — The power law fits to the angular cross-correlation 
function of lO^'^ — lO^'^ Mq galaxies with galaxies in the 10^ — 10® 
(red), 10^-10'' (green), 10^-10* (blue), 10* - 10^ M© (cyan) and 
IQf - IQiO M0 (magenta) range. 

auto-correlation functions. Figure [7] shows the cross- 
correlations of low and high mass galaxies for the Dl, 
D2 and D4 fields. Again the errors are dominated by 
the field to field variance, with the effective volume fur- 
ther reduced by the overlap integral of the two redshift 
distributions, Fig.[TJ Figure [8] shows the results of power- 
law fits to the measured angular cross-correlations over 
our redshift range. The Limber inversion to a physical 
cross-correlation, (xo/r)'', uses the redshift distributions 
of both sets of galaxies in the cross-correlation. The re- 
sulting [a;o,7] pairs are shown in Figure (5] Comparison 
to Figure [6] shows that the [xq, 7] pairs are similar to the 
[ro,7] pairs as a function of mass. The two lowest mass 
bins have switched positions with the cross-correlation of 
the lowest mass galaxies being stronger than the auto- 
correlation, although the difference is at about 2 stan- 
dard deviation level in either xq or Xg, which are the im- 
portant quantities for pair counts. The difference arises 
from the best fit, since the point-by-point differences are 
all statistically equal for the two lowest mass correlations. 

4. LOW MASS GALAXY NUMBERS 

To measure the number of low mass satellite galaxies 
around a massive galaxy we need to determine the mean 
volume density of dwarfs. Our z = 0.1 — 0.18 sample is 
entirely composed of low mass dwarfs with high star for- 
mation rates-otherwise they would be too faint. Direct 
evidence that we are missing low star formation rates 
galaxies is shown in Figure fTOl where we display the dis- 
tribution of the specific SFR in the 10^ - 10^ Mq bin, 
comparing the z = - 0.03, - 0.05 and 0.1 - 0.18 red- 
shift ranges and showing the increasing incompleteness. 
At redshift 0.02 our sample reaches Mi ~ —8 mag and 
our 0.75" resolution corresponds to 20Qh^^ pc, sufficient 
to reach ab out 100 of the 120 dwarf galaxies in the lo- 
cal volume (iBelokurov et al. 20071 iGilmore et al. 2007t 
ISharina et al. 2008f ). The simplest interpretation is that 
these three specific SFR (sSFR) distributions are con- 
sistent with a similar underlying sSFR distribution with 
increasing bias towards high specific star formation rates. 
There is a likelihood of some evolution in the sSFR dis- 
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Fig. 9. — The Limber's equation inversion of the angular cross- 
correlation to give the 3D cross-correlation length of 10^ -10^'^ Mq 
galaxies with galaxies in the 10^ — 10^ (red), lO'' — 10^ (green), 
10^-10** (blue), 10* -10^ M0 (cyan) and lO^'-lO^" Mq (magenta) 
range. 
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Fig. 10. — The volume distribution of specific star formation rate 
for the ^ = 0.1 - 0.18 (red), z = - 0.05 (green), and 2 = 0- 0.03 
(blue) samples. The SFR is averaged over 0.5 Gyr and the method 
sets a limit of about 10~*yr~^ on the maximum value. 



tribution even over this small redshift range so it would 
be unwise to assume an invariant sSFR distribution to 
estimate the total numbers of low mass galaxies. On the 
other hand, these star formation rates are not so high 
as to lead to significant evolution of the overall num- 
bers of galaxies in this mass range, so we use the den- 
sity of low mass galaxies at low redshift as our best es- 
timate. All three densities are shown in Figure [11] If 
we shrink the bin to z = — 0.02 we find the densities 
at all masses rise, which suggests that the volume is be- 
ing underestimated, about as expected given the median 
redshift error of 0.02. The range we use z = — 0.03 
does not significantly suffer from this over- estimation at 
higher mass so we believe that the densities are essen- 
tially correct. The steepening of the luminosity function 
at low mass has been remar ked for the Virgo cluster and 
other nearby systems Q iTrentham fc Hodgkin 20021 : 
ISabatini et al. 2003tlRi^es fc Geller 2008f) . 
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Fig. 11. — The field volume density of galaxies as a function 
of mass for the z = 0.1 — 0.18 (red), z = — 0.05 (green), and 
2 = — 0.03 (blue) samples. 

As a point of interest, a very simple approach to esti- 
mating the incompleteness is to use the specific SFR as 
an estimate of the duty cycle for visibility at high star 
formation rate. This will not work for high mass galaxies 
since the entire population is visible, independent of the 
star formation duty cycle. Normalizing to the highest 
mass bin, Mq, we estimate the numbers at lower mass 
as n'{M) = n{M)sSFR{M)lsSFR{MQ). The result is 
shown as the dashed line in Figure [TT] and provides a 
reasonable estimate of the numbers, at about the factor 
of two level of precision. We do expect that a duty cycle 
model should be applicable at these low masses where 
the relatively short bursts of star formation cannot be 
sustained and the sSFR distribution appears to have a 
small, fairly narrow peak at the high sSFR, with most 
of the galaxies distributed over much lower sSFR. The 
analysis of Lee et al. (2008) finds a duty cycle of 6% 
for a sample that is about a factor of ten more luminous 
than ours. Using the simple analysis above we find very 
similar duty cycles of 2.7% for 10^ - 10^ Mq and 4.1% for 
10^ — 10^ Mq. We note that the ratio of the star forming 
dwarfs to those which include the quiescent galaxies is 
6.3%, 6.8%, 4.5% in the three lowest mass bins, remark- 
ably similar to Lee et al. (2008). Again this boosts our 
confidence that the redshifts are sufficiently accurate to 
place the galaxies into the appropriate bins. 

5. FIELD GALAXY SATELLITE COUNTS 

N-body simulations (iSpringel et al. 20081 ; 

iDiemand Kuhlen fc Madau 2007t iDiemand et al. 2008( 1 
provide quantitative predictions for the numbers of 
low mass sub-halos around a Milky Way type galaxy. 
There is some disagreement over predicted numbers, 
which is likely largely due to the details of the input 
cosmological parameters, the density perturbation 
spectrum slope and normalization and the specific 
details of small numbers of realizations. Using our 
cross-correlation and density statistics we calculate 
N{M,r) = 47rn(M) ^'^{xq{M) / x)'^ x^dx, where r is set 
to the estimated 50 times critical density radius, 433kpc 
(for h=0.73) which is the radius of comparison for the 
Aquarius and Via Lactea simulations. For comparison 
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Fig. 12. — The average number of galaxies as a function of mass 
inside a radius of 433 kpc (H=73) for an adopted M/M« = 30, 
approximately M/Ly — 100, for the z = 0.1 — 0.18 (red), z = 
— 0.05 (green), and z = — 0.03 (blue) samples. We consider the 
2 = — 0.03 sample to be largely complete to M ~ —8, comparable 
to local group samples. An absolute luminosity scale for quiescent 
galaxies with Mt/Ly = 3 is shown on the top axis. The dashed 
lines are the Aquarius simulation results and Aquarius divided by 
3.1 which approximate the Via Lactea result (dot-dashed line). 



to the simulations we need to estimate the dark matter 
sub-halo mass associated with our dwarf galaxies. 
We use MjM^ — 30, which for quiescent galaxies is 
MjLy ~ 100. The counts with this normalization are 
displayed in Figure 1121 Although very low luminosity 
dwarfs can have MjLy approachi ng 1000, these ar e too 
faint to enter into our sample (jSimon fc Geha 2 007). 
Increasing the dark matter ratio a factor of three, to 
M/M* = 100, or a quiescent M/Ly ~ 300, would 
alleviate the discrepancy of Fig. [12] slightly but not 
nearly enough to make any qualitative improvement. 
Similarly if we ascribe all our dwarf s with M» < 10'' Mq 
to be resident within lO'^M© halos (jStrigari et al. 2008f) 
the discrepancy is hardly changed. 

The average number of satellite galaxies with M* > 
10^ Mq per massive galaxies is 78 and 39, for our field 
density estimates to z = 0.03 and 0.05, respectively. The 
systematic errors and field to field variance dominate the 
numbers so we will take our result to be 60 ± 20 satel- 
lites. To the same limit, which we take to be My — —8 
mag, the Milky- Way has about 30 satellites. The duty 
cycle approach gives a total of about 175 satellites on 
the average. Our assessment is that given the various 
uncertainties the agreement between our field estimates 
and the Milky Way counts is surprisingly good, with our 
estimates tending to be about a factor of two higher. 
The trend of an increasing discrepancy to lower mass is 
also seen. At a halo mass of 10^ Mq the numbers are 
a factor of 30 to 100 below the simulation results. The 
problem is one of the slope of the observed relation is 
much shallower, approximately dN/dlnM = —0.5, than 
the predicted one —0.9, not a case of some mass scale 
at which the relationships diverge. The purely random 
statistical errors of about 20% are about the size of the 
symbols. 
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6. DISCUSSION 

This paper uses faint galaxy photometric redshifts to 
find very low stellar mass galaxies that are nearby. We 
estimate that at least 80% and possibly 90% of the red- 
shifts are within the stated error of 0.02 in redshift, with 
the rest being fairly uniformly distributed over the entire 
redshift range searched, z = — 3. The analysis in this 
paper uses redshift bins that are larger than the errors 
and angular correlation analysis to minimize the impact 
of redshift errors. From The tests and comparisons in 
this paper establish the utility of the sample, but there 
is little doubt that obtaining spectroscopic redshifts and 
other follow-up studies would be of great interest. 

The correlation length for 10^ — 10^ Mq galaxies is 
about OAh~^ Mpc. There is no clear prediction for com- 
parison to our measurements. Linear biasing theory cal- 
ibrated at higher mass predicts a value about 76% of 
the linear regime, about 3/i^^ Mpc, far above our mea- 
surement. The limited evidence from existing n-body 
simulations suggests that the nonlinear bias is likely to 
be much smaller than linear theory gives, and may even 
fall into the range of our measurements. 

To estimate the number of satellite galaxies within the 
50 times critical over-density radius Milky- Way, approx- 
imately 400 kpc, requires us to make a measurement 
of the field volume density of our dwarf galaxies. The 
cross-correlation sample at a mean redshift of about 0.15 
is completely missing the low sSFR dwarfs, but they 
are visible in lower redshift samples which we there- 
fore use to estimate the density. Combining the den- 
sities and the cross-correlations predicts the number of 
low mass galaxies, to 10^ M© inside a large halo as ap- 
proximately 60 ± 20, as compared to 30 for the Milky- 
Way to a comparable limiting stellar mass. The radial 
distribution of our dwarfs is much more centrally con- 
centrated than the n-body simulation predictions which 
find that most of the objects are at large radius, whereas 
our cross-correlation. Fig. [71 is equivalent to an ap- 
proximately linear rise with radius, and the slope of 
the relation does not vary with radius. Although the 
statistics do not allow a comparison of the radial slope, 
both the Milky- Way and M31 dwarfs are also much 
more centrally concentrated than the simulated sub-halos 
dKravtsov ct al. 2004; Willman et al. 2004). 

One of the possible explanations for the deficiency of 
dwarfs is that those low mass halos that did not collapse 
before reionizat ion had their star fo rmation largely sup- 
pressed forever ([Bullock et al. 2001 1 ). The "illuminated" 
halos have numbers similar to the dwarf galaxy counts 
(|Strigari et al. 2007f l. The mass of dwarf galaxies is ev- 
erywhere dominated by dark matter with the stars con- 
tributing relatively little to the binding energy. Stellar 
astrophysics adds energy to the gas which normally works 
to help drive gas out of the dwarf. The dwarf galaxy 
stellar populations are a m ix of stars of every age, with 
significant ongoing bursts (jMateo 1998f ). The evidence 
suggests that virtually all dwarfs cycle through bursting 
and passive stages (Lee et al. 2008). By extension, it 
would not be possible to stop low mass dark halos that 
did not happen to have stars formed near reionization to 
form stars at a later time. An observational argument is 
that the suppressed star formation would lead to a range 
of M/Ly values that extended to infinity, whereas in our 



range of absolute magnitude more luminous than -8 mag 
the distribution appears to be well bounded. Hence, the 
census of low mass dark halos is likely fairly complete 
within our mass limit. 

It is possible that the total extent of the dark matter 
is much larger than the st ellar contents of a dwarf galaxy 
U iPenarrubia et al. 20081 ) in which case the mass will be 
underestimated. A dark matter mass to luminous mass 
ratio of approximately 1000, roughly M/Ly — 3000, 
at a stellar mass of 10^'^ Mq is required to agree with 
the Aquarius simulation, decreasing to dark to luminous 
mass ratio of 30 at a stellar mass of 10*'^ Mq. Such 
very large M/Ly values have not been advocated in the 
observational literature, although the mass outside the 
stellar distribution could lead to such large values. Our 
field density, about 5 Mpc"'^, of dark-to-stellar mass ratio 
1000 halos imply an flM ^ [0.01,0.03,0.03,0.10] in the 
bins with mean mass [lO^ "*, lO^ ^, lO^ "*, 10*-^] Mq which 
comes close to the total dark matter density but does not 
exceed it. 

The fact that the radial distribution of dwarfs are more 
centrally concentrated than the simulations imposes a 
dynamical limit on the M/Ly of the satellite galaxies. 
Either the numbers need to rise by about 30 or the 
masses need to go up by about 30 times to bring the to- 
tal counts inside the virial radius into basic agreement. 
However, both of these (if uncorrelated) will increase 
disk heating by a factor of 30, leading to considerably 
more disk heating than the shallo wer distributions find 
(jFont et al. 200lUArdi et al. 2003f ) likely at a level that 
would be unacceptable. 

Our overall conclusion is that we have found approxi- 
mately twice as many dwarf galaxies per galaxy, 60 ± 20, 
as the 30 in the same stellar mass range as for the Milky- 
Way. For conventional M/Ly ~ 100 (a dark-to-stellar 
mass ratio of 30) this then slightly alleviates the missing 
satellite problem, but still leaves the numbers about a 
factor of 30 short of the n-body predictions. Although 
there are no clear predictions for the two-point auto- 
correlation function in the non-linear regime for low mass 
galaxies, the very short correlation length of OAh~^ Mpc 
is smaller than expected from the linear regime. The 
outcome of this paper is to extend the measurements of 
dwarf clustering to the field which results in a strength- 
ening of the discrepancy with LCDM clustering expecta- 
tions. 

Peebles (1989, 2007) has pointed out that dwarf galax- 
ies, as a low bias population, should be a reasonable 
tracer of the dark matter and should therefore be present 
in the voids between massive galaxies, whereas it appears 
that the dwarf galaxies largely follow the massive galax- 
ies. The fields appear remarkably uniformly populated 
with dwarf galaxies. Figure [3l however it is important 
to recognize that we only measure clustering to approxi- 
mately OAh~^ Mpc, not the 10/i~^Mpc scale and larger 
that characterize voids. 
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